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Abstract 

In this study, we investigate the reflectance property of the cylinder, right circular cone, and square pyramid 
shapes of silicon nitride (Si3N4) subwavelength structure (SWS) with respect to different designing parameters. In 
terms of three critical factors, the reflectance for physical characteristics of wavelength dependence, the reflected 
power density for real power reflection applied on solar cell, and the normalized reflectance (reflected power 
density/incident power density) for real reflectance applied on solar cell, a full three-dimensional finite element 
simulation is performed and discussed for the aforementioned three morphologies. The result of this study shows 
that the pyramid shape of SWS possesses the best reflectance property in the optical region from 400 to 1000 nm 
which is useful for silicon solar cell applications. 



1. Introduction 

Silicon solar cell is one of the promising renewable 
energy technologies in order to relieve the impact of the 
climate change. In semiconductor-based solar cells, elec- 
tron-hole pairs are generated through absorption of 
impinging photons. Due to high refraction index of 
semiconductor materials, especially silicon, the incident 
sunlight power is largely reflected back, resulting in the 
reduction of light absorption and poor energy conver- 
sion efficiency. Antireflection coating (ARC) is mounted 
over absorption layers, resulting in three effects: (a) 
reduction in surface reflection, (b) increase in light 
absorption due to an increase in optical path length by 
diffraction, and (c) enhancement of internal reflection 
that reduces the amount of escaping light. Based on the 
theory of impedance matching, single layer (SLR) and 
multilayer of ARC are proposed for reduced reflectance 
property; however, the resulting reflectance spectra meet 
the demand only within a narrow spectral domain. Sub- 
wavelength structure's (SWS) dimensions are much 
smaller than the wavelengths of light; therefore, using 
ARC on the surface of silicon solar cells can substan- 
tially reduce the reflectivity and improve the capability 
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of light trapping. It thus will achieve the enhanced effi- 
ciency according to our recent both numerical and 
experimental studies [1-3]. Compared with silicon solar 
cell with a SLAR, the efficiency of silicon solar cell with 
Si3N4 SWS is promising among various ARC layers in 
our recent work [4]. A rigorous coupled-wave analysis 
(RCWA) [1,5-7] has been reported to estimate the 
reflectance of Si3N4 SWS by approximating structural 
shapes with partitioned uniform homogeneous layers. 
RCWA is an exact solution of Maxwell's equations for 
the electromagnetic diffraction by grating structures 
which is generally applicable for 2D plane with ID peri- 
odicity; however, RCWA may suffer numerical difficul- 
ties in presence of evanescent orders and it requires a 
large amount of calculation for retaining several diffrac- 
tion orders. These factors limit flexible application of 
RCWA; in particular, for 3D problems with non-azi- 
muthally symmetric structural shapes. Numerical simu- 
lation of 3D morphological effect on reflectance 
property has not been studied yet. Therefore, a full 3D 
finite-element (FE) analysis of Si3N4 SWS will be an 
interesting examination for quantitative understanding 
of the reflectance property. 

In this study, 3D FE simulation for the reflectance of 
Si3N4 SWS with three types of structural shapes, the 
cylinder, the right circular cone, and the square pyramid 
shapes, is conducted with respect to different geometry 
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parameters and lighting angles for quantitative under- 
standing of reflectance property. First, proper selection 
on the boundary conditions can alleviate the computa- 
tional load from simulating a holistic ARC. The reflec- 
tance of Si3N4 SWS on the silicon substrate is thus 
simulated using the 3D finite element method (FEM); 
consequently, in terms of three critical factors, the 
reflectance for physical characteristics of wavelength 
dependence, the reflected power density for real power 



reflection applied on solar cell, and the normalized 
reflectance (reflected power density/incident power den- 
sity) for real reflectance applied on solar cell are calcu- 
lated and discussed for the aforementioned three 
morphologies. The analysis of reflectance spectrum with 
wide-angle incidences of electromagnetic wave and the 
average reflectance with various heights are presented. 
Besides, according to our recent study, which presented 
the optimal design parameters of Si3N4 SWS based on 
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Figure 1 (a) Plot of the periodic structure of Si3N4 SWS with 1x1 and 2x2 arrays as unit cell. 3D schematic plots of the examined (b) 
cylinder-, (c) right circular cone-, and (d) square pyramid-shaped structure, respectively. 
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RCWA [4], numerical verification and comparison is 
accomplished following the discussion. The engineering 
findings of this study show that the pyramid shape of 
SWS possesses the best reflectance property in the opti- 
cal region from 400 to 1000 nm which is useful for sili- 
con solar cell applications. 

This rest of the article is organized as follows. In Sec- 
tion 2, we show the computational structure and model. 
In Section 3, we report the results and discussion. 
Finally, we draw conclusions and suggest future work. 

2. The SWS and optical model 

Based upon our experimental characterization, Figure la 
illustrates a periodical structure of Si3N4 SWS which is 
used in our 3D FE simulation without loss of generality. 
We study Si3N4 SWS with the cylinder, the right circu- 
lar cone, and the square pyramid shapes, as shown in 
Figure Ib-d, respectively. With a constant volume, the 
diameter of cylinder- and right circular cone-shaped 
Si3N4 SWS and the edge length of square pyramid are 
130 nm, the heights (h) of the etched part of Si3N4 SWS 
are 200, 600, and 471.3 nm, the height (s) of the non- 
etched part is 70 nm, and the base (W) of a unit cell is 
200 nm [4]. Note that the thickness of Si substrate is 
given 600 nm. Note that all structural parameters are 
adopted from our experimental studies [2-4,8]. 
Throughout the article, we consider time-harmonic 
fields assuming a time-dependence in e'^^\ The diffrac- 
tion problem is governed by the well-known Maxwell 
equations 



E2 (f)=£i(f + L)=£i (f)e- 



(5) 



V X E = —jcofjiH, 

V X H = J + jcosE, 

V -E = p/s, 
and 

V • B = 0, 



(1) 



(2) 



(3) 



(4) 



where E and D are electric field intensity and flux 
density, H and B are magnetic field intensity and flux 
density, A is the corresponding frequency to the wave- 
length A, /, and p are current density and charge density, 
s is electric permittivity, fi is magnetic permeability. A 
repeated pattern is applicable to use periodic boundary 
conditions, thus the Floquet theorem is adopted to 
simulate the boundary condition of periodic structure. 
Floquet theorem asserts that the analysis region can be 
reduced significantly in one periodicity cell to character- 
ize the propagation property. The electric fields in peri- 
odic structure are related as follows: 



where r is position vector, L is the distance between 
the periodic boundaries, and ^ is a phase factor deter- 
mined by wave vector k and L: 



0 = k-L, 



(6) 



The polarization of transverse electric (TE) mode, in 
which the electric field is normal to the direction of 
wave propagation, is excited as the normal incident light 
source with wavelengths sweeping from 400 to 1000 
nm. The bottom region of Si substrate is assigned as 
perfect matched layer in avoidance of reflected wave. 
The refraction index of Si3N4 is 2.05, and the refraction 
index of Si is frequency dependent with the relation [1]: 
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where A is the incident wavelength, A ■■ 
8.10461 X 10 ^ A.1 = 1.1071 (im, and s 



0.939816, B = 
■■ 11.6858. The 
calculation settings of reflectance were reported and can 
be found in our recent studies [1,4]. 

3. Results and discussion 

In order to examine the effect of Floquet boundary con- 
dition in 3D FE analysis, as shown in Figure 2, we com- 
pare the difference between the simulated unit cells of 1 
X 1 and 2x2 array of Si3N4 SWS. We find at the wave- 
lengths above 600 nm, the reflectance of 1 x 1 array of 
Si3N4 SWS as unit cell is almost consistent with unit 
cell of 2 X 2 array, meanwhile insignificant discrepancy 
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Figure 2 Plot of the difference of reflectance spectrum of Si3N4 
SWS with the cylinder-, the right circular cone-, and the square 
pyramid-shaped structures as well as two different periodical 
configurations: 1x1 (solid line) and 2x2 arrays (dashed line) 
in the 3D FEM simulation. 
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Figure 3 Comparison of tiie reflectance spectra for the cone- 
and cylinder-shaped Si3N4 SWS calculated by RCWA and 3D FE 
analysis with the same design parameters 



occurs at wavelengths shorter than 600 nm. Based on 
this consequence, it is enable us to do simulation with 
more computational efficient using 1x1 array of Si3N4 
SWS as a simulated unit cell with engineering accepta- 
ble accuracy. According to our recent RCWA work [4], 
the reflectance spectra are first plotted in Figure 3 using 
the optimal design parameters [1,4]. Also, the spectra 
calculated by a full 3D FE analysis with the same design 
parameters are indicated by dashed lines. For the cylin- 
der-shaped Si3N4 SWS, the reflectance spectra for 
RCWA and FE analysis are similar, but not agreed for 
the cone-shaped Si3N4 SWS due to existing evanescent 
orders along the top of structures. This comparison 
confirms the importance of 3D FEM simulation which 
is beyond the RCWA approach [4]. 

Figure 4a-c shows the reflectance spectra with incident 
angles of 0°, 15°, 30°, 45°, and 60° for the cylinder-, right 
circular cone-, and square pyramid-shaped Si3N4 SWS, 
respectively. For the normal incidence case, the lowest 
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Figure 4 Plots of the reflectance spectrum for the (a) cylinder- (b) circular-cone-, (c) and square-pyramid-shaped Si3N4 SWS with 
incident angles of 0°, 15°, 30°, 45°, and 60°. 
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average reflectance among three structural shapes is 
3.47% of square pyramid-shaped structure. The others 
are 6.86 and 4.42% for the cyUnder- and the right cir- 
cular cone-shaped Si3N4 SWS, respectively. Meanwhile, 
as shown in Figure 4a-c, one can observe that the 
reflectance increases significantly with larger incident 
angles, resulting in average reflectance beyond 50%. 
Table 1 summarizes the average reflectance for various 
incident angles. Height effect on average reflectance of 
Si3N4 SWS at normal incident angle with d = 130 nm 
and 5 = 70 nm is also calculated, as shown in Figure 5. 
The resulting average reflectance of pyramid-shaped 
Si3N4 SWS nearly keeps lowest in comparison with the 
cylinder- and the right circular cone-shaped Si3N4 
SWS as the structural height is ranging from 50 to 500 
nm. Figure 6 shows the reflectance dependence on the 
structural height and wavelength. The pyramid-shaped 
Si3N4 SWS has lower reflectance and less sensitivity on 
structure height in comparison with the cylinder- 
shaped Si3N4 SWS. Hence, the impact of process varia- 
tion of structure height on solar cell performance is 
smaller for pyramid-shaped Si3N4 SWS. Based on solar 
spectrum at the sea level revealed in American Society 
for Testing and Materials (ASTM) Standard Tables for 
Reference Solar Spectral Irradiances: Direct Normal 
and Hemispherical 37 Tilted Surface [9], we further 
estimate the reflected power density (W/m^/nm) 
defined by reflectance times incident power density, as 
shown in Figure 7. The higher reflected power density 
of cylinder-shaped Si3N4 SWS (red line) indicates the 
less efficiency in the solar cell application. Therefore, 
normalized reflectance defined as 

_ Reflected power density 

^^norm = " 1 7— (8) 

Incident power density 

reveals the real power efficiency applied in the solar 
cell application. Figure 8 shows the normalized reflec- 
tance for the cylinder-, right circular cone-, and square 
pyramid-shaped Si3N4 SWS, respectively. The square 
pyramid-shaped Si3N4 SWS again shows the lowest nor- 
malized reflectance 3.13% while the cylinder- and the 
right circular cone-shaped Si3N4 SWSs have 6.66 and 
4.12%, respectively. 



Table 1 Summary of the average reflectance of SI3N4 
SWS with various incident angles 
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Figure 5 Plot of the average reflectance among the studied 
three shapes of Si3N4 SWS with heights varying from 50 to 
500 nm. 




400 500 600 700 800 900 1000 

Wavelength (nm) 



(a) 




400 500 600 700 800 900 1000 

Wavelength (nm) 



(b) 

Figure 6 3D view for the height effect on the reflectance with 
respect to different wavelength, (a) The pyramid-shaped Si3N4 
SWS has lower reflectance and less sensitivity on structure height in 
comparison with (b) the cylinder-shaped Si3N4 SWS. 
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Figure 7 Plot of the reflected power density among three 
different shapes. 



4. Conclusions 

In this study, the reflective property of unit cell with a 
validated Floquet boundary condition has been calcu- 
lated using a full 3D FE simulation. Considering various 
incidence angles and height effect on three experimen- 
tally observed structural shapes of Si3N4 SWS, we have 
concluded that the pyramid-shaped Si3N4 SWS has best 
reflective property in the analysis of morphological 
effect. Compared with the results of RCWA, the reflec- 
tive property calculated by the full 3D FEM is signifi- 
cantly deviated from the results from RCWA, giving the 
hint that a detailed and comprehensive methodology is 
dispensable for the design of Si3N4 SWS. The results of 
computed reflectance, reflected power density, and nor- 
malized reflectance have shown that the pyramid shape 
of SWS may have the best reflectance property in the 
optical region from 400 to 1000 nm which is useful for 
silicon solar cell applications. The optimized pyramid- 
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Figure 8 Plot of reflectance with and without considering 
incident solar spectrum at sea level. 



shaped Si3N4 SWS is currently under plan for imple- 
mentation with silicon solar cells. 
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